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Abstract 



We present a theory of particles, obeying intermediate statistics 
Q ' ("anyons"), interpolating between Bosons and Fermions, based on the princi- 

ple of Detailed Balance. It is demonstrated that the scattering probabilities of 
identical particles can be expressed in terms of the basic numbers, which arise 
, naturally and logically in this theory. A transcendental equation determining 

^ ' the distribution function of anyons is obtained in terms of the statistics pa- 

' rameter, whose limiting values and 1 correspond to Bosons and Fermions 
00 ' 

. respectively. The distribution function is determined as a power series involv- 

OO , ing the Boltzmann factor and the statistics parameter and we also express the 

\ distribution function as an infinite continued fraction. The last form enables 

C3 ' one to develop approximate forms for the distribution function, with the first 
approximant agreeing with our earlier investigation. 
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I. INTRODUCTION 



We formulate a theory of particles obeying intermediate statistics, interpolating between 
Bosons and Fermions, which might be called anyons. Our formulation will be based on two 
assumptions: 1) The exchange symmetry or permutation of the coordinates of the particles 
in the many particle wave function is accompanied by multiplication by a complex number 
/, thus generalizing the symmetric or antisymmetric wave functions and 2) The principle of 
Detailed Balance: if ni,n2 represent the average occupation numbers of states labelled by 
1 and 2, then the number of transitions flowing from 1 to 2 must equal that flowing from 2 
to 1 at equilibrium. 

The particles described by this theory may or may not be the same as the anyons obeying 
intermediate or fractional statistics discussed in the hterature. The objects named anyons 
[1-8] carry both an electric charge and a magnetic flux. They have attracted a great deal of 
attention and have been the subject of intense investigation in the literature. The anyons 
arise from the special circumstance in 2+1 space-time dimensions, where the permutation 
group is the braid group and the Chern-Simons theory provides a natural realization of the 
anyons [9,10]. There has been a great deal of discussion in the literature on the thermo- 
statistics of anyons [11]. Since the real world is strictly in a 3+1 dimensional space, anyons 
may not be real particles: they could be quasi-particles playing important roles in condensed 
matter phenomena. More recently, the subject of generalized statistics has been investigated 
in one dimension in the context of many kinds of statistics [12]. In contrast to the theory of 
anyons familiar in the literature, our present approach is not limited to two space dimensions 
and is valid in the real world of 3+1 dimensional space-time. It is important to point out 
that this is due to the fact that we do not invoke the spin-statistics theorem and do not 
require ordinary spin to be interpreted in two dimensional space. 

On one hand, the theory based only on Detailed Balance does not have the features of the 
full-fledged theory of anyons which takes advantage of the braid group in two dimensions. 
On the other hand, in this formulation we investigate the idea of interpolating statistics only 
in the context of the statistical mechanics of a gas in equilibrium, without the constraints 
imposed by quantum field theory. Here we shall use the name anyons, for convenience, to 
refer to particles obeying intermediate statistics, with the understanding that the present 
work needs to be developed further by incorporating additional assumptions before the 
connection with true anyons described in the literature could be established. 

The subject of anyons has been well investigated in the literature, especially in the 
context of quantum field theory and the braid group [2,5,7,8]. Interesting results have 
also been derived to describe the thermostatistics of anyons, such as determining the virial 
coefficients [11,13]. However, the theoretical basis of the statistical mechanics of anyons has 
not really been established. For instance, the distribution function for the anyons has not 
been determined in an exact form. We had used an ansatz for the distribution function 
in our earlier work in order to derive many of the thermodynamic properties [11]. More 
recently, Chaturvedi and Srinivasan [14] have done a comparative study of the different 
interpolations between Bose-Einstein (BE) and Fermi-Dirac (FD) statistics available in the 
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literature, in the context of the most general interpolations, including Haldane statistics and 
Gentile statistics, with the conclusion that the distribution function introduced by us [11] 
has some desirable features. 

Consequently it is worthwhile to investigate the basic theoretical structure of anyons 
from the point of view of statistical mechanics and investigate the distribution function 
of the anyons if possible. This is the goal of this paper. Not only shall we determine 
the distribution function of the particles obeying the interpolation statistics, we shall also 
formulate a theory which leads to this determination without resorting to the restriction to 
a two dimensional space. Our formulation will be based on the quantum theory of many 
particles permitting a generalized interpolating exchange symmetry statistics with no other 
assumptions. Remarkably, this formulation leads to an exact theory which requires the 
employment of basic numbers. 

In Sec. II we study the scattering probability of many particle states obeying interpolating 
statistics. We establish the fact that the basic numbers arise naturally and logically in this 
theory. We introduce the method of Detailed Balance in Sec. Ill and derive a transcendental 
equation for the distribution function in terms of the Boltzmann factor and the statistics 
determining parameter. The distribution function of anyons is studied in Sec. IV where we 
show that a closed form solution is not possible. We present the exact solution for the dis- 
tribution function as an infinite series as well as in the form of an infinite continued fraction 
which is amenable to approximations. Sec. V contains a brief summary and conclusions. 

Let us introduce an ensemble of particles, the anyons, which obey a generalized statistics, 
interpolating between (BE) and (FD) statistics. We begin with the framework for building 
the wave functions of these anyons by a generalized procedure of /-symmetrizing in such a 
way that it will reduce, in appropriate limits, to the standard procedure of symmetrizing for 
Bosons and anti-symmetrizing for Fermions. The operation of pcrnmtation or exchange of 
the coordinates of the many particle wave function results in multiplication by the complex 
number /, the exchange symmetry parameter, so that 

^*n(- • • , • • • , • • •) = /*n(- ■■,qj,---,qi,---)- (1) 

Since the Hamiltonian has the property V~^7iV = 7i, it follows that TiiV'^n) = EniV'^n) 
and thus V^n is an eigenfunction of the Hamiltonian with the same eigenvalue as Thus 
f^n is proportional to ^„ and consequently we may take |/P^n — Hence the general 
exchange symmetry that would lead to an intermediate statistics could be implemented by 
the complex number / with the property |/p = 1. We shall accordingly choose the exchange 
symmetry to be implemented by / = e*'^", where a is the statistics determining parameter, 
< a < 1, so that /* = f"^. The hmits / — > 1,-1 correspond respectively to the 
BE and FD statistics representing Bosons and Fermions. This procedure for incorporating 
exchange symmetry among anyons is justified [7] by the special property of rotations in two 
dimensions. Here we may treat it as an assumption or an ansatz, not restricted to two space 
dimensions but valid in any number of dimensions. 
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II. MANY PARTICLES AND QUANTUM PROBABILITIES 



Following Feynman [15], we consider the two particle scattering amplitude defined by 
the product 0162, where ai = (l|a) describes the scattering of particle a into state 1 and 
62 = (2 1 6) describes the process 6^2. We shall take 1 and 2 to be the same state at the 
end in order to deal with identical particles. The exchange symmetry has to do with the 
process corresponding to a — 2, 6 — 1 which is indistinguishable from the direct process 
and the amplitude for this process would be / a2 bi due to the exchange factor /. The 
total probability amplitude is the sum of the direct and exchange processes. Employing 
the abbreviation (l|a) = (2|a) = a, we find the probability of this two particle scattering 
process involving non-identical particles to be 

pii=(i+i/ni«n&r = (2) 

since |/p = 1. This probability is the same as for ordinary Bosons [15]. However, for 
identical particles, we need to take account of interference between the two processes and 
that makes a great deal of difference. We obtain the probability in this case to be 

p2Lw = (i+i/r+/+r')i«n&r, (3) 

since /* = and this probability depends on the statistics determining parameter a. In 
the limit / — > 1, this would reduce to the case of Bosons and it would be twice as much 
as in Eq.(2) for the non-identical particles. For arbitrary / the probability for the process 
involving identical particles relative to that for non-identical particles is given by 



p(2) 
identical 



i(2 + / + /-^) \am. (4) 



In what follows, we shall omit |ap, |6p etc. since the single particle states are properly 
normalized. Understanding this to be a proportionality, we might henceforward refer to this 
itself as the probability for the process, or just the probability of the two particle state, for 
simplicity. 

Similarly, we consider the three-particle processes a— >1,6— >2,c— >3 together with 
the exchange processes with factor / for each exchange operation, thus resulting in the 
combination abc + / act + / bac + bca + cab + cba with the same abbreviation as 
earlier. We can now determine the probability for the three particle process, or just for the 
three particle state with identical particles, as 

p(3) ^1(1 + 2f-' + 2/-2 + /-=^)(1 + 2/ + 2f + f) . (5) 

We can reduce this to the form 

p(3) = 1 |io + 8(/ + /-I) + 4(f + + if + r')} 

= ^ {6 + 7(/ + r') + 4(f + r' + 1) + if +r'+f+ r') } ■ (6) 
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Recognizing the pattern here, we observe that the right hand side in the expressions for the 
probabilities in Eqs.(4), (5), (6) contain basic numbers [16] , with the base /. They are 
indeed expressed succinctly in terms of the basic numbers defined by 

We shall henceforward omit the subscript / for simplicity. Here / = e^^°' , the BE hmit 
is q; — >^ 0, / — s> 1 and the FD limit is a — 1,/ — —1 and /* = f~^. Our formulation is 
symmetric under / and is the familiar generalization of the basic numbers introduced 

long ago by F. H. Jackson [16]. Studying the limits we find that the Bose limit gives [n] — > n 
while the Fermi limit is quite different: [n] ( — 1)""''^ n which becomes —n for even numbers 
and +n for odd numbers. In this limit, we are therefore dealing with a generalization of 
ordinary fermions for which the exclusion principle is not automatically valid: in the FD 
limit it can be imposed by hand only. Returning to the basic numbers introduced in Eq.(7), 
it is quite evident in Eqs.(4), (5), (6), and as will be seen in more detail below, that the 
basic numbers arise naturally, automatically, in the theory of particles obeying interpolating 
statistics. This feature is quite new and not recognized in the standard theory of anyons 
in the literature. This might be an indication that the theory of interpolating statistics 
naturally involves a deformation of the system such as that described by the basic number 
system with its consequences, without however introducing a deformed algebra of operators. 

In addition to the representation, 

[n] = r-' + r-' +■■■ + /-"+^ + /-"+^ , (8) 

the following representation of the bracket number 

r , sinnTTo; 

W = , 9 

smTTo; 

will be found most useful. In terms of the basic numbers, the probabilities may now be 
expressed conveniently and succinctly as 

pP) = i(2 + [2]), 

p(3)^ l(6 + 7[2]+4[3] + [4]) , 

p(4) ^ 1 ( 35 + 54[2] + 52[3] + 36[4] + 18[5] + 6[6] + [7] ) , (10) 

and so on. We can derive a useful result, 

[l] + [3] + [5] + --- + [2n-l] = ([n]f . (11) 
which can be proved by using [17] the identity. 



1 'A . .^ /sinnt\^ 



5 



We shall also find the following result quite useful: 

[n - 1] [n] = [2] + [4] + [6] + ■ ■ ■ [2(n - 1)] . (13) 

In other words, [2] [3] = [2] + [4], [3] [4] = [2] + [4] + [6], etc. This is proved by using the 
identity 

• , . {n-l)y . ny y 

} ^ smky = sm sm — esc-. (14) 

k=o 2 2 2 

These identities can be used to re-express Eqs.(lO) in the form 

p(^) = i(2 + [2]), 

p(3) = i.(2 + [2])(2 + 2[2] + [3]), 

P(^) = i (2 + [2]) (2 + 2[2] + [3]) (2 + 2[2] + 2[3] + [4]) . (15) 

In this manner, generalizing to n particles, we are led to the probability for the n-anyon 
state: 

p(") = ^ (2 + [2]) (2 + 2[2] + [3]) • • • (2 + 2[2] + 2[3] + • • • + 2[n - 1] + [n]) . (16) 
n ! 



III. DETAILED BALANCE 

From Eq.(16), we can infer the enhancement factor, which is a measure of how much 
greater the probability of the n + 1 -particle state is, compared to the probability of the 
n-particle state. For now, we simply express it as 

p(n+l) 

nn) = , (17) 

which can be determined from Eq.(16). This enhancement factor provides the essential step 
in the method of Detailed Balance. This brings us to the second assumption upon which our 
theory of interpolating statistics rests: if ni,n2 represent the average occupation numbers 
of states 1 and 2 respectively, then the number of transitions flowing from 1 to 2 must equal 
that flowing from 2 to 1 at equilibrium. This is the principle of Detailed Balance. We 
should stress that this principle is characteristic of thermodynamic equilibrium and may be 
regarded as a consequence of the second law of thermodynamics [18]. Indeed the principle 
of Detailed Balance is valid when thermodynamic equilibrium prevails or the validity of 
microscopic reversibility in the language of statistical physics. This notion is based on 
the reversibility of the microscopic equations of motion, or on the Hermitian nature of the 
scattering Hamiltonian [19]. The principle of Detailed Balance can thus be stated as 

niF(n2)e^^i = n2F(ni)e^^2 ^ (18) 
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where the population of each level is governed by the Boltzmann factor and F{n) is the 
enhancement factor. This yields 



e^-^ — constant — z , (19) 

F{n) 

which would, in principle, enable us to determine the distribution function for the anyons 
in terms of the Boltzmann factor. In the BE limit for instance, the enhancement factor is 
just n + 1, which immediately leads to the BE distribution 

n = — - . (20) 

where z = e^^ is the fugacity of the gas. Let us now proceed in the same manner for 
arbitrary /. The enhancement factor for arbitrary / is given by Eqs.(16) and (17) which 
reduces to the simple form 

= = ;^(2 + 2[2] + 2[3] + • • • 2[n] + [n + 1]) . (21) 

To proceed further, we employ the result 

n 

^ [A;] = 2 cos(7ra/2) [n/2] [(n + l)/2] , (22) 

fe=0 

which is easy to prove by using identities involving sums of trigonometric functions. This 
gives us the following result for F(n), after some algebra: 

F{n) = ^— - 2 cos(7ra/2) [(n + l)/2] { [n/2] + [n/2 + 1] } . (23) 

We can also derive another identity 

[n/2] + [n/2 + 1] = 2 cos(7ra/2)[(n + l)/2] , (24) 
which can be used to simplify the above result. Hence we obtain 

^^""^ = ^^^"^ + ^^/^^ cos(7ra/2)}^ . (25) 

We observe that this reproduces the expected result F{n) — > n + 1 in the Bose hmit. Upon 
now invoking the Detailed Balance we obtain the important result 

-/"^ = F{n)/n = ^^^^^{[{n + l)/2] cos(W2)}^ ■ (26) 

This can be rewritten in the following convenient form in order to deal with the distribution 
function: 

g/3(E-/.) ^ 1 sin=^(n + l)7ra/2 

n(n + l) sin^7rQ;/2 ' ^ ^ 



Solving this equation should, in principle, lead to the distribution function for the anyons 
in this formulation. At this point before dealing with the distribution function, we need to 
study the nature of the intermediate statistics as an interpolation between the BE and FD 
hmits. Specifically we need to study the hmits corresponding to BE and FD statistics. 

We have already observed that F{n) ^ n + 1 in the Bose limit, / ^ 1. Indeed, it is 
readily verified that Eq.(26) reproduces the correct BE statistical distribution. The case of 
Fermi limit, a— >!,/—>— 1, is, however, somewhat complicated. We have: [n] —>■ n,—n 
for odd and even occupation numbers respectively. Evaluating the limit in Eq.(26), the 
enhancement factor may be put in the form, 

lim F(n) = ^— (im e^("+i)'^/2)' , (28) 

thus resulting in a generalized Fermion theory. This corresponds to an infinite dimensional 
representation analogous to the generalized fermions investigated by Chaichian et al . If we 
consider only n = 0, 1, and impose the exclusion principle "by hand", then F{n) 1, 0. This 
is equivalent to F{n) — 1 — n and in this case it reproduces the standard Fermi distribution 
with the exclusion principle. For arbitrary n, which is allowed, however, it is an infinite 
value representation 

limF(n) = ^{l,0,l,0,l,---}, (29) 

and we find that F(n) has repeating values l/(n + 1) and in the limit a — * 1, / — — 1. 
This is thus a special and interesting feature of this theory. Thus for arbitrary values of 
the parameter /, including the Fermi limit / — > — 1, the theory of intermediate statistics 
requires the existence of generalized fermions beyond the exclusion principle. 



IV. THE DISTRIBUTION FUNCTION 

We begin by rewriting Eq.(27) in the form 

g,(^-.) ^ 1 sin>+l)x _ 
n{n + 1) sin x 

Here x — T:a/2 in terms of the statistics determining parameter. The object is to determine 
the average occu})ation number n, the distribution function, so that we can compare with the 
standard BE or FD distribution and understand the nature of the interpolating statistics. 
We can expand the right hand side in a power series 

e'^^^~^''i = i + ao + ain + a2n^ + asn^ H , (31) 

n 

where 
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ao = csc^ x{x sin 2x — sin^ x) 

Qi — csc^ x{x'^ cos 2x — X sin 2x + sin^ x) 

02 = csc^ X {~^^ ^'^^ ~ ^^^^ X + (x — 2x^/3) sin 2x| 

03 = csc^ X — x^/3) cos 2x + sin^ x — (x — 2x^/3) sin 2x| , (32) 

and so on. The coefficients to any desired order can he ol)tained by using Mathematica. It 
is clear that ao — > 1, while a„ — > for n > 1 in the Bose limit, which is consistent with 
the BE distribution as described in Sec. III. In the case of the Fermi hmit, Oq — > —1 and it 
leads to the generahzed fermions. We can rewrite the above as 

1 ai 2 "^2 3 as 4 /oo\ 
- — n n n n — ■ ■ ■ . [66) 

9 9 9 9 

where g = e^^^~^^ — oq. This series can be reverted to express n as a series in powers of 
thus 

n^l/g + B/g' + C/g' + D/g' + E/g'' + F// • • • , (34) 

where 

B = ai/g, C = {2al/g^ + 02/^), D = {5aia2/g^ + a^/g + 5al/g^) , 
E = E = eaiflg// + 802// + 14ai// + (^4/ 9 + 2X0102// , 

F = 70104/5^ + 70203/51^ + 840^02// + 05/51 + 2801O3// + 2801O2// + 42oi// , (35) 

and so on. Upon rearrangement, we can rewrite the above form for the distribution function 
as the following series: 

n{E) = 1/g + ai/g^ + 02// + (2o? + 03)// 

+ (5oi02 + 04) // + (5o? + 601O3 + 3o^ + 05)// H . (36) 

We observe that this is a power series in 1/gr but the term with l/g"^ is absent. We can 
rewrite Eq.(36) in the form 

n{E)^l/g + f2ak/g', (37) 

fc=3 

where 

<^3 — cii, 0(4 — a2, 2ai + 03 = 0:5, 5oi02 + 04 = 0:6, (38) 

etc. We can now express this in the form of a continued fraction. We invoke the standard 
algorithm which can be introduced as follows [21]. If a continued fraction is of the form 

Bj, 61 / „s 

T^-Co + 4 , (39) 

C2 H 
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then the successive convergents (approximants) are obtained by the recurrence formulae 



Bn = CnBn-l + 6„i?„_2, C„ = C„C„_i + 6.„C„_2, = 1, C_i = . (40) 

Thus in order to put the series in Eq.(37) in the continued fraction form, we proceed with 
the first convergent and set cq = Bq/Cq = 0. Using the recurrence relations, Eq.(40), 
we determine Bq = 0, Co = = l,ci = g,Bi = l,Ci = g. This determines the first 
convergent as n = 1/g. Next, the recurrence relations lead to 62 = —0(3g, C2 — g^ + 0^3, B2 — 
g'^ + as, C2 = g^ which determines the second convergent to be 

n = . (41) 

9 - — 

5' +0:3 

Continuing on with successive convergents in this manner, we obtain the desired infinite 
continued fraction form for the distribution function as follows 

n{E) = (42) 

9 

2 , oi4g 
g 



oi^g 

a^g + 0^4- 



a^g + as 

This is the distribution function in the exact theory, albeit not in a closed form. Approxi- 
mations can be implemented in order to investigate the thermostatistical properties of the 
anyons. The first approximant, 

n(E) = } , (43) 

was indeed employed, as an ansatz, in our earlier investigation [11]. It might be pointed out 
that the continued fraction form is mathematically desirable in view of the fact that it can 
be expressed as a ratio of polynomials, hence analytic in the variable a as well as offering 
the best approximation possible. This form is also amenable to investigation in terms of 
Pade approximants. 

Finally it must be pointed out that Eq.(42) is an exact expression for the distribution 
function, albeit in the form of an infinite continued fraction, and can be determined explicitly 
for any specific value of the parameter a, either as an infinite continued fraction, or to any 
desired approximation. For instance, for the case of a = |, the exact form of the distribution 
function is 

n{E)Ui = , (44) 

^ ooS^ 

^ 2 ^ {ao-nyA8)g 
g -ao + 



aog - (ao - 7r748) + 



where ao(|) = n/2 — 1. 
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V. SUMMARY 



The thermodynamic distribution function of anyons in two space dimensions, obeying in- 
terpolating statistics is not known in the hterature and is an open question. We have not only 
found an answer, thus determining the distribution function, we have also demonstrated that 
the theory describing interpolating statistics has several remarkable and interesting features. 
Our investigation leads to a generalized definition of permutation symmetry in arbitrary di- 
mensions and not restricted to two space dimensions. We have shown that the theory of 
permutation symmetry that would describe particles obeying interpolating statistics is suc- 
cinctly formulated in the language of basic numbers. These basic numbers arise naturally 
and automatically in this formulation but do not explicitly invoke any deformed oscillator 
algebra. Our theory is based on the principle of Detailed Balancing which is a consequence 
of the Second law of thermodynamics. Furthermore, this formulation leads to the determi- 
nation of the exact distribution function, without having to introduce any approximation. 
In this manner, we have formulated the theory which leads to a transcendental equation for 
the distribution function of anyons in terms of the statistics determining parameter and the 
Boltzman factor containing the energy of the state. We obtain a solution of this equation 
which we express as a power series as well as in the form of a continued fraction. We show 
that the first approximation of this theory reproduces a form of the distribution function 
introduced by us in an earlier investigation. 

An important feature of our formulation consists of the fact that the basic numbers arise 
naturally and automatically in this theory, specifically the symmetric formulation of the 
basic numbers. The basic numbers arise automatically in our theory but they are not part 
of an oscillator algebra, nor do we introduce the construction of Fock space for the particles 
obeying the intermediate statistics. The theory reduces to BE statistics in the Bose limit, 
/ ^ 1,« ^ 0. The Fermi limit, / —l,a 1 leads to generalized fermions beyond 
the exclusion principle and thus correspond to infinite dimensional representation and they 
reduce to the familiar FD statistics only when n is restricted to 0, 1 by hand. 

Although the algebra of q-oscillators appears in the literature on the subject of anyons 
[7,8,20], it is also known that q-oscillators may have nothing to do with anyons since the for- 
mer exist in arbitrary space-time dimensions. In our formulation, we do not use q-oscillator 
algebra, neither do we use Fock states. We do not make use of the Hamiltonian other than 
to recognize that the Hamiltonian should be of a form that permits Detailed Balance. The 
Hamiltonian consists of only the kinetic energy terms and is the free Hamiltonian [5] in the 
anyon gauge and the anyon wave functions satisfy "twisted" boundary conditions. 

It is interesting that the first approximant (first convergent) of our solution corresponds 
to the approximate form for the distribution function 

n{E) = / — , (45) 

which is the form introduced in an earlier investigation [11]. It is interesting that this also 
agrees with the partition function derived by the method of Green function [20] , to first order 
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approximation. It should be stressed that while the form above is a point of agreement for an 
approximate theory, the form of a in the present formulation as in Eq.(32), as a function of 
the statistics determining parameter, is very different from that of [11]. Finally we obtain the 
exact form of the distribution function for the case a = | in the form of a continued fraction 
which can be evaluated to any desired order. Such a form for the distribution function can 
indeed be obtained for any chosen value of a, the statistics determining parameter 

It is of much interest that this formulation leads to an intermediate statistical mechanics, 
interpolating between BE and FD, without explicit use of the special properties valid in 
2+1 space-time dimensions, but only stems from the principle of Detailed Balance. This 
contrasts with investigations in the literature of anyons which arise from the braid group in 
the special case of 2+1 space-time. As is seen from Eq.(30), the dependence on the statistics 
determining parameter is through a quadratic function, hence invariant under a — > —a, 
implying a clockwise-counterclockwise symmetry of the braid [6] . It would be worthwhile to 
continue this investigation further in order to explore the relation to these anyons. As this 
formulation reveals several desirable and interesting theoretical features, it is appropriate to 
pose some questions such as the following. Does the generalization to any space dimensions 
imply that the connection with Chern-Simons type of gauge theory can be extended beyond 
two space dimensions? Is it possible to formulate the partition problem corresponding to 
this generalized, exact theory? Is it possible to determine the many useful thermodynamic 
properties of the system described by such an exact theory? These questions will occupy us 
in further investigations of this theory. 
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